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Abstract: In this paper, we firstly consider the properties of the Shannon entropy for two probability 
rS^ [ distributions which stand in the relationship of majorization. Then we give a generalized Uhlmann theo- 

rem in an infinite dimension Hilbert space. Also, we show that S(<&(p)) — S(p) for all quantum states p if 
and only if there exists an isometry operator V such that $(p) = VpV* , where $ is a quantum channel. 
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1 Introduction 

m 

Let B(H) be the von Neumann algebra of all bounded linear operators on a separable Hilbert space 
H and S(H) be the set of all density operators on %. That is p £ S(H) if and only if p > and 
tr(p) = 1. Each element of S(H) is also called a quantum state in quantum information. Let T(H) be 
the set of all trace class operators on %. An operator V £ B(H, /C) is called an isometry if V*V = /%. 
In this case, VV* £ B(JC) is an orthogonal projection. As usual, for x,y £ H, x <S) y denotes the 
operator z — > {z,y)x (z £ %). For p £ S(H), we denote \(p) = (Xi(p), ^(p), • • • A n (p) • • •), which 
is the sequence of eigenvalues of p with the non-increasing order. Thus X(p) £ Cq, where Cq is the 
positive cone of sequences decreasing monotonically to as in [7]. Let Z°°(M) and /{(IR + ) denote the 
set of all bounded real sequences and all summable non-negative real sequences, which have sum 1, 
respectively. For a vector r £ l°°(M), let H- = (r|, . . . r„ • • •) denote the vector whose elements are the 
elements of r re-ordered into non-increasing order. Adopting the definition of majorization in [2,7], 
for r,s £ Cq, we say that r is majorized by s, written as r -< s, if 

k k oo oo 

XM < XM' for fc = l,2,--- and^r| = ^s}- 

8=1 j=l i=\ i=l 
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Recently, many mathematicians have paid much attention to the infinite dimensional Schur-Horn 
theorem and infinite majorization. In [11], A. Neumann has given properties of infinite majorization 
in Z°°(R) and V. Kaftal and G. Weiss in [7] have obtained interesting results for infinite majorization 
in Cq. Arveson and Kadison by different methods presented some other characterizations in [3]. Some 
related results are the following: if r, s G /J(M + ), then 

r -< s -4=> r = Qs, with Qij = \Uij\ for some unitary U [6, Theorem 1], 

and 

r -< s -4=>- r = Qs, for some orthostochastic matrix Q [7, Corollary 6.1]. 

Motivated by the above results, we firstly consider the properties of the Shannon entropy for two 
elements in Z^(1R + ) which have majorization. Then we extend and study those properties for two 
operators p,a G S(7i). Following the definition given for finite dimensional spaces (see [1]), we denote 
p -< a for two operators p, a G S(H) if X(p) -< X(a). 

Let M n {B{%)) be the von Neumann algebra of n x n matrices whose entries are in B{T-L) and 
$ : B(H) — > B{U) be a linear map. Then $ induces a map $ n :M n (B{H)) — ► M n (B(H)) by the 
formula 

$„(Ki)) = (Q(aij)), for (aij) G M n (B(H)). 

If every $ n is a positive map, then $ is called completely positive. $ is said to be normal if $ is 
continuous with respect to the ultraweak (a— weak) topology. Normal completely positive contractive 
maps on B{T~L) were characterized in a theorem of Kraus [8, Theorem 3.3], which says that $ is a 
normal completely positive map if and only if there exists a sequence {ylj}^ 1 of B{%) such that for 
all X G B{H), 

oo oo 

<f>(X) = ^TAiXA* with ^AiA*<I, 

i=l i=\ 

where the sequence {Aj}"^ is not necessarily unique and X)i=i A i A t — I m the strong operator 
topology. The family {Aj}?^ is also called a family of Kraus operators for $. In this case, the dual 
of <I> is defined by 

oo 

&(X) = J2 A *i XA i forXGT(^). 
i=i 

It is easy to see that \tr[&(X)Y]\ = \tr[X$(Y)]\ < \\$\\\\Y\\tr(\X\), for X G T{H) and Y G B{H), 
so &(X) G T{T~L) is well defined on T{%). In general, 3>t cannot be an extension from T{T~L) into 
B(H). However, if X^i A*Aj < I, then 3>t is well defined on B{T~L) and <J?t i s normal. A normal 
completely positive map <3? which is trace preserving ("&'(/) = I) is called a quantum channel. If a 
normal completely positive map satisfy <!>(/) < I, then <I> is called a quantum operation ([8,9]). A 
quantum operation is unital if &(I) = /, which is equivalent to ^ • AjA* = I. A quantum operation 
is bi-stochastic if it is both trace-preserving and unital. In particular, $ is said to be a mixed 



unitary operation if &(X) = ^2™ = itiUiXU* , where n < oo, Ui are all unitary operators and £« > 0, 

The von Neumann entropy of a quantum state p is defined by the formula 

S(p) = -tr(plog(p)) 

As usual, in this formula logarithms are taken to base two. In classical information theory, the 
Shannon entropy is defined by H(p) = — ^2 i Pilog(pi), where p = (pi,P2 • • • Pn ■ ■ ■) is a probability 
distribution. If A» are the eigenvalues of p, then the von Neumann entropy can be re-expressed as 

oo 

i=i 
where we use log = 0. 

In [4], Hardy, Littlewood and Pdlya showed for £,77 G M. n , 

£ -< f] <^=^ £ = Qf] for some doubly stochastic matrix Q. 

In the quantum context, Uhlmann obtained 

p -< a -4=/- p = $(cr), $ is a mixed unitary quantum operation, 

for p,a G S(7i), where % is a finite dimension Hilbert space. This Uhlmann theorem might be used 
to study the connection between majorization and quantum mechanics. Furthermore, the quantum 
operations which can preserve the von Neumann entropy and the relative entropy of quantum states 
were extensively studied in more recent papers [5,10,15,17]. 

In this paper, we firstly consider the properties of Shannon entropy for two probability distribu- 
tions which have majorization. Then we give a generalization of the result due to Uhlmann [1] in 
an infinite dimensional Hilbert space. Also, we show that S($(p)) = S(p) for all quantum states p 
if and only if there exists an isometric operator V such that 3>(p) = VpV* , where <3? is a quantum 
channel. 

2 Shannon entropy of infinite probability distributions 

The following lemma is a direct corollary of [13, Theorem 8.0.1] as the function f(x) := — xlog(x) is 
a concave function. 

Lemma 2.1. Let a, b G R n with a», h > 0, and a -< b. Then Ya=i fi a i) > Ya=i f( b i)- 

The following is an extension of above lemma from the finite case to the infinite majorization. 

Theorem 2.2. Let a, b G l\(R + ) and a ■< b. Then H(a) > H(b). 



Proof. Suppose that c = (ci, c 2 , ■ ■ ■ Cn, • • ■), where q > and Xl£i Q < 1. It is clear that for 
arbitrary n, we have 

n 

(Ci, C2, • • • Cn, 0, 0, 0, • • •) -< (£ Q, 0, 0, 0, • • •), (2.1) 

1=1 

so Lemma 2.1 implies — Y17=i c * 1°§ c « ^ ~(Y17=i c «) l°g(Ej=i c j)- Letting n — > 00, by the continuity 
of the function f(x) = — xlogx, we get 

00 00 00 

- y, ^ iogQ > - y, c * lo sE c j)- ( 2 - 2 ) 

i=l j=l j=l 

Firstly, we assume that b has only finite nonzero elements. Without loss of generality, we suppose 

(ai, a 2 , ■ • • a n , ■ ■ ■) = a -< b = (h, b 2 , ■ ■ ■ , b m , 0, • • •), (2.3) 

where bi > 0, for 1 < i < m. Thus there exists N such that X/i=iV+i a * < min{6i, 62, ■ ■ ■ &m} ; so 

00 
(ai,a2,---ajv, ^ a i> >°>---) "< (bi,62,---,^m,0,0---), (2.4) 

i=JV+l 

which implies 

N 00 00 m 

-^ajlogaj- ^ Ojlog( ^ clj) > -^6jlog(6j). (2.5) 

i=l 1=^+1 j=JV+l i=l 

By inequality (2.2), we know H(a) > H(b). 

For the general case of 6, we note that for any positive integers s, the relation a -< b entails 

00 
(ai,a 2 ,---a n ,---) -< (bi,b 2 ,--- ,b s , Y ^AO---), ( 2 -6) 

i=s+l 

so the proof above yields 

S OO OO 

H(a) > - Y bdog h - Y bi log( £ h). (2.7) 

1=1 j=S+l j=S+l 

Letting s — > 00, we get H{a) > H(b), as desired. 

Lemma 2.3. Let a,b G R n with c^A > 0. If a ■< b and ELi/K) = E*=i Z( 6 *) where 
/(x) = — xlog(x), then a^ = b^. 

Proof. Without loss of generality, we assume that 

a 1 = (01 , 02 , ■ • • a m , • • • 0) ■< It = (61, 62, • • • , b m , • • • 0), (2.8) 

where m < n and aj + 6j 7^ 0, for 1 < i < m. For < t < 1, denote 

c t = (ten + (1 -t)bi,ta 2 + (1 - t)b 2 ,---ta m + (1 - £)&m)- 



It is easy to verify that a -< ct -< b, so Lemma 2.1 implies 

m 

- Y,( ta i + 0-- Wi) l °s{tai + (1 - t)bi) = H(a). (2.9) 

i=i 

Taking the second derivation function at t of the both sides of the equation(2.9), we get 

f-' tai + (1 - t)6i 
«=i 

which yields Oj = 6j, for 1 < i < m, so or = hr. 

The following proposition shows that the Shannon entropy of the infinite probability distribution 
is strictly monotone in the relation of majorization. 

Proposition 2.4. Let a, b G ^j(K + ) such that all a^ > 0. If a -< 6 and H(a) = #(6) < oo, then 
a* = b^. 

Proof. Let 

a 1 = (ai,a 2 , ■ ■ ■ a n , ■ ■ ■) -< 6+ = (61,62, • • • , 6 n , ■ • •)■ 

Then a\ < 61, so there exists /c such that ai S (6fc+i,6fc], which implies a\ = tb^+i + (1 — £)&&, for 
some < t < 1. Denote 



c = (61, 6 2 , ■ ■ ■ %+i + (1 - *) fc fc) C 1 ~ & fe+i + tb k , b k +2, ■ ■ ■), 

it is clear that a -< c -< b. Thus by Theorem 2.2 and the assumption, we have H(a) = H(b) = H(c), 
so Lemma 2.1 implies a\ = tbk+i + (1 — t)bk = bk- Let 

a = (a 2 , a 3 , ■ ■ ■ a n , • • •) and 6 = (61, 6 2 , ■ ■ ■ 6 fc _i, b k+1 , ■■■). 

Thus a -< b and -ff(a) = -ff(6) < 00, so by a similar proof, we conclude that there exists k\ > k such 
that a 2 = 6fej . By mathematical induction, we get a n+ \ = bk n , for some subsequence of b. However 
YaL\ a i = SSi bi = 1 and all aA > imply a n+1 = b kn = b n +i, that is a 1 = 6-K 

Remark 2.5. The conditions of all Oj6j > and -ff(a) < 00 are essential in Proposition 2.4. 
Indeed, it is obvious that both conditions a ~< b and H(a) = H{b) < 00 are not changed if we add 
some zeros for a. Furthermore, if H(a) = H(b) = 00, we may replace a by a' = (^-, -?f, a 2 , ■ ■ •)• Then 
a' -< b and H(a!) = H(b) = 00. However, it is a contradiction that ar = b^ = a!^. 

3 von Neumann entropy of quantum states 

Let us recall that two quantum states p\ and p 2 are /^-equivalent if and only if there is a sequence 

of unitary operators {Ui}^ 1 such that lim„ >00 ||pi — U n p 2 U*\\i = 0. 

Proposition 3.1. Let pi,p 2 G S(H). If p\ -< p 2 and S(p\) = S(p 2 ) < 00, then pi and p 2 are 
C 1 — equivalent. 



Proof. By Proposition 2.4, all the non-zero spectral points (including the multiplicity) of p\ and 
P2 are the same, which is equivalent to the fact that p\ and p2 are /^-equivalent, by [3, Proposition 
3.1]. 

The following two results give extensions of Uhlmann theorem in an infinite dimensional Hilbert 
space. We need some eigenvalue estimates of Weyl given in [16]. Let A be a compact operator 
with eigenvalues Ai > A2 > A3 > • • • and V n be the set of all n-dimensional projections. Then 
max P67 ,„ fr(AP) = X)?=i A*. 

Proposition 3.2. Let pi,p2 G S(7i). If p\ is finite rank, then p\ -< p2 if and only if there exists 
a mixed unitary operation <£ such that <&{p2) = Pi- 

Proof. Necessity. If p\ is finite rank, then denote A(pi) = (Ai(pi), A2(pi) • • • A m (pi), 0, • • •). As 
Pi -< P2, we have that p2 is also finite rank and A(p 2 ) = (Ai(p 2 ), A 2 (p 2 ) ■ ■ ■ A n (p 2 ), 0, • • •) where n < 
m. By the spectral decomposition of states p± and P2, we conclude that there exist two orthonormal 
bases {x{\ c *l 1 and {yj}£i OI " % such that 



Pi = ^ \{pi)xi <g> Xi and p 2 = ^ \{p2)yi <8> y%- 



i=l i=l 



Let Tii be the subspace spanned by {xi}™ =1 and U be the unitary operator defined by Uyi = Xi, for 
i = 1, 2 ■ ■ ■ . Then pi|^j and Up2U*\-n 1 can De represented bymxm matrices, which have spectra 
(Ai(pi), A 2 (pi) • • • A m (pi)) and (Ai(p 2 ), A 2 (p 2 ) • • • A n (p 2 ), • • • 0) G M m , respectively. Thus Uhlmann's 
theorem ([1 or 13 Theorem 4.1.1]) implies that there exists a mixed unitary operation <I> on S('H 1 ), 
that is Pil-Hj = YLlLi UVi(Up2U*) \-HiV*, where Vi G B(H 1 ) are unitary operators, XX=i ^* = 1 an< ^ 
all tj > 0. As (J7p2i7*)|^x = 0, we denote Ui = diag(Vi,I^±), so l7j G H('H) are unitary operators, 

and pi = £^^([^2 £/*)£/*. 

Sufficiency. By the Weyl estimates of [16], we get for n = 1, 2 ■ ■ ■ , 

E"=i A i(Pi) = m axp e p„{Tr(piP)} 

= maxpep^TrCEZiWiPtUiP)} 

= ma,x PeVn {Y,T=i t i Tr (p2U*PU i )} 

< Z?=iU™x PePn {Tr(p 2 U*PU t )} 
= max P6 p n {Tr(p 2 P} 

= YJl=l^i{p2)- 

Theorem 3.3. Let pi,p2 G S(H). Then the following three conditions are equivalent: 

(a) pi -< p 2 ; 

(b) there exists a sequence of mixed unitary operation ^> n and a bi-stochastic quantum operation 
* on S(H) such that linin^oo ||* n (p) - ^(p)||i = for all p G S(H), and *(p 2 ) = pi; 

(c) There exists a bi-stochastic quantum operation ^ such that ^(p2) = Pi- 

Proof, (a) ==£- (6). If pi -< p2, then A(pi)^ -< A(p 2 ) , it follows from [6, Theorem 1] that 
A(pi)^ = QX(p2)^, where Q is an infinite matrix satisfied Qij = \Uij\ 2 for a unitary operator U. 
By the spectral decomposition of states P2, there exists an orthonormal basis {yjj^i of % such 



that p 2 = Yl'iLi ^i{P"i)Vi ® Vi- Let {xj}^ 1 be the usual orthonormal basis of H (in the sense of 
isomorphism), that is X{ = (0, • ■ ■ 0, 1, ■ ■ •). Denote U\ = Y^=i x * ® £/« an d e « = U^U*Uiyi. Then it 

V v ' 

i 

is easy to see that {ei}'^2 =1 is an orthonormal basis of % such that Xi(pi) = (p2&i, &i) for i = 1, 2 ■ ■ ■ . 
Thus p2 has an infinite matrix form as (with respect to the basis {ej}?^) 



( Ai(pi) A12 
A21 A 2 (pi) 



P2 



Vil 



A In 
A2n 

A n (pi 



... \ 



V 



(3.10) 



• / 



where Ay = Ajj. Denote the sequence of rank-one projections Ei = a <8> e^, then define 3>(p) = 
Si=i EipEi, for p G S(H). Also for n = 1, 2 • • • , let $ n be the mixed unitary operation defined as 

follows <J>„(p) = \ Ya=i UipiU*)*, where J7 = diag(u,uj 2 ■ ■ ■ u 11 , 1, 1 ■ ■ ■ 1 ■ ■ ■) and u> = exp *n . By 
a direct calculation, we get 

1 n ( e, (8) ej, i = j 

I 0, i 7^ j and i < n or/ < n. 



n 



i=l 



Thus 



where 



$n{p2 



P21 
P22 



:(\/W)©( V < e *»> 

i=n+l 



i=l 



(3.11) 



(3.12) 



P21 = dia5(Ai(/?i),A 2 (pi)---A n (pi)), /> 22 = (/ - P„)p 2 U - -Pn)|(V£ n+1 M), 
and P n denotes the orthogonal projection on the subspace Vr=i{ e *i- Then 

||$ n (p 2 ) - $(p 2 )||i = 2Tr[$ n ( / 9 2 ) - $(/> 2 )] + < 2Trp 22 — >• 0, (n — ► 00), 
since Tr& n {p2) = Tr&(p 2 ) implies 

||$n(p 2 )-$(P2)||l = Tr(|<l> n (p2)-<J>(P2)|) 

= Tr[$ n (/9 2 ) - $(/> 2 )]+ + Tr[$ n (p 2 ) - $(p 2 )]" 
= 2Tr[$ n ( P2 )-<J>(p2)] + , 

and $ n (p 2 ) - $(p 2 ) < (/ - ^n)p 2 (/ - P„) yields [$ n (p 2 ) - $(/> 2 )] + < P+(/ - PJfri! ~ P n )P+, so 

Tr[$ n (p 2 ) - $(p 2 )}+ < Tr[P+(I - P n )p 2 {I - P n )P+] 

< Tr[(I - P n )p 2 {I - P n )] 
= Trp 22 , 

where A + , A~ are the positive and negative parts of the self-adjoint operator A, and P+ is an 
orthogonal projection on the range of [Q n {p2) — ^{p2)} + • By the spectral decomposition of p±, we 



conclude that there exists an orthonormal basis {/i}^i of T~L such that p\ = YH^=i ^i(pi)fi ® /«• 
Define a unitary operator V by V/j = ej, for i = 1,2---, so pi = F*<&(p 2 )V. Denote ^n(p) = 
y*$ n (p)F and *(p) = F*<I>(p)V, for all p G £(?{); then by the proof above, we similarly get 
lim„^ 00 ||* n (p) - *(p)||i = lim n _^oo ||$n(p) - $(p)l|i = 0, for all p G S(H). (b) =► (c) is clear. 

(c) => (a). Let pi = \l/(p 2 ) = Y^S=\ AiP2A*. Applying Weyl's eigenvalue theorem, we get for 
n = l,2---, 

n oo oo 

J>(pi) = max{Tr(^^p 2 4*P)} = max {Tr(p 2 J>*PA)}- (3.13) 

As £*=i A*A i = I and £^ A A* = J > then E*=i A* p 4 < J and 



Tr£>*PA) = Pr(f>A*^) 



n. 



For convenience, we denote B = ^2i^=i-A-*PAi and Ui = (Bxi,Xi), where {xi}°2 =1 is an orthonormal 
basis of H satisfying p 2 = ESi ^i{P2)xi <& X{. Thus < U{ < 1, and E£i Uj = n, so 

Tr(p 2 B) = 52™ i Ai(p 2 )ui = £? =1 Ai(p 2 )ni + £~ n+1 A;(p 2 )ui 

< EJLi MP2)w; + A„+i(p 2 ) E£n+i u * (3 u \ 

< E?=l A i(P2)nj + A n+ i(p 2 ) YJt=i( l ~ U i) 

< TJUHP2). 

By equation (3.13) and inequality (3.14), we conclude that for all n = 1, 2---, J27=i^i(pi) — 
X)"=i ^i(P2) as desired. 

Remark 3.4. In a finite dimensional Hilbert space, p\ -< p 2 is equivalent to p\ = 3>(p2), for some 
mixed unitary operation $. However, for an infinite dimensional Hilbert space, p\ -< p 2 does not 
imply p\ = ^2i?-i tiUip2U* , where Yl'hLi ^i = lj U > an d Ui are unitary operators for all i. Indeed, 
the condition p\ = Y^Li UUifcU* yields dim(ker(pi)) < dim(/cer(p 2 )). But we can supplement many 
zeros for A(pi). 

The following proposition was obtained in [10,17] for the finite case, in which the condition of 
injectivity of p may be dropped. 

Proposition 3.5. Let p G 5("H) and $ be a bi-stochastic quantum operation. If p is injective 
and S(p) = S(Q(p)) < oo, then 3>(p) = UpU* for a unitary operator U. 

Proof. Suppose 3>(p) = Xli^i^iMt*- Then by Theorem 3.3, we have 3>(p) -< p. We claim that 
if p is injective, then so is ^(p). We assume that $(p) were not injective, then there is a vector 
i /0 that satisfies <3?(p)x = 0, so (X^i AipA*x, x) = 0, which yields A*x = 0, for all i. Thus 
x = E^i AiA*x = 0, which is a contradiction. It follows from Proposition 2.4 that A(p)^ = A($(p))^, 
so by the spectral decomposition theorem, we get ^(p) = UpU* for a unitary operator U. 

Remark 3.6. In an infinite dimensional Hilbert space, the condition that p is injective may not 
be dropped. 

In the following, we shall characterize the structure of a quantum channel that does not change 
the von Neumann entropy of any quantum state. In [11], Molnar and Szokol gave the structure of the 



map which preserves the relative entropy in a finite-dimensional Hilbert space. However, we cannot 
find the paper which study the structure of the map which preserves the von Neumann entropy. The 
following two lemmas are needed. As usual, /C('H) is the set of all compact operators on T~L. 

Lemma 3.7. (see [14]) (i) For all p G S(H), S(p) > and S(p) = if and only if p is a pure 
state. 

(ii) In a n dimension Hilbert space, S(p) < logn, and S(p) = logn if and only if p = —. 

Lemma 3.8. (see [9]) Let $ be a quantum operation with &(I) < I. Then {B G K,(H) : ®(B) = 
B} = {B G K,{U) : &{B) = B} C A', where A' is the set of commutators of A = {A^A* : A L G 
B(H),i = 1,2- ■■}. 

Theorem 3.9. Let $ be a quantum channel on B{T-L). Then S(Q(p)) = S(p) for all quantum 
states p G S(H) if and only if there exists an isometry V G B(H, JC) such that ®(X) = VXV* for all 
X €B(H). 

Proof. Sufficiency is clear. 

Necessity. For any unit vector x, we have S($>(x x)) = S(x (£> x) = 0, so <l?(x x) is a rank one 
orthogonal projection. Thus there exists a unit vector z such that <&(x (8) x) = z z. 

Let y-Lx be another vector of %. For convenience, denote 



<fr(x ® x) = x' <8) x' and <J>(y y) = y' <g) y'. 



Setting po = x 9 y , we have 



g ^®^'®^ ) = 5($(p )) = S(p ) = I- (3.15) 

In the following, we shall show x' _L y' . Let %o ^ "H denote the two-dimensional space spanned by x' 
and y'. Then 3>(po) can be treated as an operator from T-Lq into %$. Suppose y' = ax' + y/l — |a| 2 x /_L , 
where < \a\ < 1 and x /_L is a unit vector in T~Lq orthogonal to x'. Then we have the following matrix 
forms for the operators x' x' and y' <S> y' , respectively: 



/„ , f 1 0\ , , ( \a\ 2 ay/1- \a\ 2 

\ U U J \ a\J\ — \a\ z 1 — \a\ 

where a is the complex conjugate of a. Then 



ay/l-\a\ 2 l-|q| 2 

2 2 



(3.16) 



By a direct calculation, we conclude that the characteristic polynomial of (3.16) is 

1 - \a\ 2 
A 2 -A + P~ = 0. (3.17) 



Furthermore, Lemma 3.7 (ii) and equation (3.15) imply that equation (3.17) has two equal roots 
Ai = A2 = i, which yields a = 0. Thus, x'_Ly' as required. From the proof above, we know that the 
map $ sends orthogonal pure states to orthogonal pure states. Let {e{\ c ?l 1 be an orthonormal basis 
of % and P n be the orthogonal projection onto the subspace spanned by {ej}™ =1 . Then 

n 

$(/) = lim $(P n ) = lim V $(e, <g> e,), 

n — >oo n — >CG ^— ' 

i=\ 

so <&(/) is an infinite dimensional orthogonal projection. 

Let p G S'(H) be injective and p = X^i ^ x « ® x *» where {xj}^ is an orthonormal basis and Aj 
are the eigenvalues of p. Then 

n n 00 

3>(p) = lim 3>(N AjXj ® Xj) = lim N Aj<3?(xj <g> xA = } AjX^x'j. 

i=l j=l i=l 

Denote V p Xi = x\ for all i, then $(p) = V^pV^. 

By the Kraus Theorem, we have $(p) = X^i ^P^*, where Aj G £(%), and ^~ 1 A*A, L = J, as 

<3? is trace-preserving. Then 

00 

1=1 

which yields 

00 

^f;Ap4*^ = p, (3.18) 

i=i 

for all injective p £ S(l-L). Furthermore, X^i=i ^*^i = ^"(-0 — ^1 so 

00 00 00 

J^ V^AiA^Vp < I and ^ A^VpV^Ai < £ A* A, = I. (3.19) 

i=l j=l i=l 

Then Lemma 3.8 and equation (3.18) imply 

00 00 

pY, v;a % aiv p = £ f;Ap4*^ = p, 

i=l i=l 

and 



00 00 

^ A- W4 = J>*v^*A = p, (3.20) 



«=i i=i 

so 



i = X>?WA, (3.21) 



i=i 
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which yields 



J2A*V p V p *A i = I = J2A*A i 



i=i 



i=i 



Thus 



i=i 

which implies A* (I - V p V*) = 0, that is 

A* = A* V P V*, for all injective p £ S(H). 
Using equations (3.18) and (3.20), we get 

J2 A*V p (J2 V p *A jP A*V p )V p *A t = p, 



(3.22) 



(3.23) 



(3.24) 



so equation (3.23) implies 



7 J A i AjpA^Aj 



h:i 



(3.25) 



which is equivalent to $T o $(p) = p, for all injective p G S(H). 

For any unit vector x, we claim that there exists a sequence of injective states p n £ S(7l) such 
that p n — >w* x ® x. Indeed, let {/i}i^i be an orthonormal basis of the orthogonal complement 
subspace of x, and p n = (1 - ±)x ®i|i Si=i(^/» ® /*)■ Tnen nm «— >oo 7>(|Pn - x <g> x|) = 0, so 
limn—** Tr[(p n - x x)X] = 0, for all X G B(W). Thus 



x<g)x = lim p n = lim <!>* o $(p n ) = <J>t o $i 



X (&) x 



which says 



^ A*^- (x ® x)A*Aj = & o $(x <g> x) = x<g>x 



- ',j 



for all rank one projections x® x. Then by Lemma 3.8 again, we have A*Aj{x <g> x) = (x <g) x)A*Aj, 
which implies that for &lli,j = 1, 2 ■ ■ ■ , ^4*Aj = Ajj/ and A*Ai = A^I, so A^ > and ^°° A^ = 1. It 
is clear that 



( A\\ 

A\ 



l* 



\ '■ J 



A, Ao ■■■ A. 



/ A\A X 
A\A X 

A* n A\ 
V ••• 
/ A n / 

A21I 



A\A 2 

A* 2 A 2 

KM 



A12/ 

A22/ 



Anil A n 9-f 



4* 4 



4* 4 



Aln-f 



... \ 



... / 
\ 
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Denote 



M 



I An 


A12 


Aln 


A21 


A22 


A2n 


X n l 


A n 2 


^nn 


V ••• 







so Xji = Xij and 



Aj,' I < XaXjj, for 1 < i, j 



since M > 0. Further equation (3.25) implies ^- • A*A{A*Aj = /, which yields 



*j J 



(3.26) 



£ I a« | 2 = 1 = £> 



'-.7 



Then by a direct calculation, we get 



/ j I A«j I — / y XaXjj, 



so equation (3.26) implies 



Thus 



for j = 1, 2 ■ ■ ■ . We denote V, 



\.- I 2 - 
A„ I - 



XaXjj, for 1 < i,j < 00. 



a*^i = x n i = yfx^x~e^- w n, 

so Vj are isometric operators for j = 1, 2 • • • , which yields 



V*V, = e^ 



so we get 



-1(9 



-ie,'T/* 



H. Then ^*ViVi* = e^- 1 ^* and VjV*V x = eS~ w iVj. Hence 
V x V?VjV*V x V? = VyVf and J^W*^* = VjV*, 

ViV?(I - VjV*) = and VjV? (I - VyV?) = 0, 



then ViV{ = VjV*, for j = 1, 2 ■ ■ ■ . Thus for all X G B(ft), 

^i * V, J v? 



which implies 



Xjjvjv? 



$(X) = A X XA\ + A 2 XA* 2 + • • • + A n XA* n + ■■■ = ViXVf. 



The following result is clear, by Theorem 3.9. 

Corollary 3.10. Let $ be a bi-stochastic quantum operation. Then S(&(p)) = S(p) for all 
quantum state p £ S(H) if and only if there exists a unitary matrix U such that <&(p) = UpU^ . 
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